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By modelling the coronal structures by "slowly" evolving 
Double-Beltrami two-fluid equilibria (created by the interac- 
tion of the magnetic and velocity fields), the conditions for 
catastrophic transformations of the original state are derived. 
It is shown that, at the transition, much of the magnetic en- 
ergy of the original state is converted to the the flow kinetic 
energy 
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I. INTRODUCTION 

The latest TRACE and SOHO/EIT observations have 
brought the solar corona into sharp focus. The obser- 
vations reveal: 1) the structures that constitute the so- 
lar corona are in constant motion; they are full of fast- 
moving gases, and are heated primarily at their foot- 
points (base) very close to the solar surface. The heating 
occurs in few minutes in the first ten to twenty thou- 
sand kilometers above the surface, i.e, in a rather small 
fraction of the bright part of the anchored structure. In 
direct contradiction to the predictions of some theories, 
the heating is neither uniform (throughout the loops) 
nor does it happen preferentially near the top. A di- 
rect quote sums up the situation aptly: "Moreover, not 
only heat is deposited low down, but the gas is often ac- 
tually thrust upward very rapidly. It does not merely 
'evaporate' into the coronal structures, it is often ac- 
tually thrown up there. Exactly how that happens is 
still a puzzle" (l), 2) the loops are composed of clus- 
ters of filamentary structures which are not, as believed 
before, static bodies supported by interior gas pressure 
and heated along their lengths |^ . They fill and drain so 
quickly that the gas in them must be moving nearly bal- 
listically (see latest TRACE news and e.g. (|l[) along the 
substructures, rather than being "quiescently heated". 
From a detailed study of the loops with different charac- 



teristic parameters one concludes that the heating pro- 
cess is quite non-uniform ||]. 

Transient brightenings, with their associated flows of 
cool and hot material, are also a very common phe- 
nomenon in the TRACE movies. These relatively fast 
(violent) happenings vary from small events in the quiet 
Sun to major flares in active regions; brightenings which 
are more than 10^ km apart often occur within the same 
exposure that typically lasts for 10 to 30 s This kind 
of a coincidence in the events at distant locations is sug- 
gestive of fast particle beams propagation along separate 
magnetic loops which come together at the flaring site. 
The flaring sites are generally assumed to be reconnection 
sites although observations have not establish a causal 
connection: "Direct evidence for reconnection in flares is 
difficult to find, despite the fact that it is thought to be 
the primary process behind flares" §. It is remarkable 
that often the post-flare loop systems begin to glow at 
the TRACE EUV wavelengths without substantial dis- 
tortion: reconnection that probably took place appears 
to be (largely) completed by the time the loops are de- 
tected. 

These observations pose a new challenge for the theo- 
ries of quiescent as well as not so quiescent coronal struc- 
tures and events. In this paper we examine the conjecture 
that the formation and primary heating of the coronal 
structures as well as the more violent events (possibly 
flares, erupting prominences and coronal mass ejections 
(CMEs)) are the expressions of different aspects of the 
same general global dynamics that operates in a given 
coronal region |Q . It is stipulated that the coronal struc- 
tures are created from the evolution and re-organization 
of a relatively cold plasma flow emerging from the sub- 
coronal region and interacting with the ambient solar 
magnetic fleld. The plasma flows, the source of both 
the particles and energy (part of which is converted to 
heat), in their interaction with the magnetic fleld, also 
become dynamic determinants of a wide variety of plasma 
states; it is likely that this interaction may be the cause of 
the immense diversity of the observed coronal structures 
1^-^. Preliminary results from this magneto-fluid ap- 
proach reproduce several of the salient observational fea- 
tures of the typical loops: the structure creation and pri- 
mary heating are simultaneous - the heating takes place 
(by the viscous dissipation of the flow kinetic energy) in 
a few minutes, is quite non-uniform, and the base of the 
hot structure is hotter than the rest. 
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We plan to extend the scope of the magneto-fluid the- 
ory beyond the creation of the semi-quiescent coronal 
structures by seeking answers to the following: a) can 
the basic framework of this model predict the possibil- 
ity of, and the pathways for the occurrence of sudden, 
eruptive, and catstrophic events (such as flares, eruptive 
prominences, CMEs) in the solar atmosphere, b) does 
the eventual fate, possibly catstrophic reorganization, of 
a given structure lie in the very conditions of its birth, 
c) is it possible to identify the range and relative val- 
ues of identifiable physical quantities that make a given 
structure prone to eruption (flaring), d) will an eruption 
be the result of the conversion of excess magnetic energy 
into heat and bulk plasma motion as is generally believed 
to happen in the solar atmosphere 0- jlj] ? 

We begin by identifying the quasi-equilibrium state 
of a typical coronal sturcture with a slowly changing 
Double-Beltrami (DB) state (one of the simplest, non- 
trivial magnetofluid equilibrium) . The slow changes may 
be due to changes in the sun which affect the local mag- 
netic fields, the interaction of various nearby structures, 
or disturbances in the solar atmosphere. The parameter 
change is assumed to be sufficiently slow that, at each 
stage, the system can find its local DB equilibrium (adi- 
abatic evolution). The slow evolution must conserve the 
dynamical invariants: the helicity hi, the generalized he- 
licity and the total (magnetic plus the fluid) energy 
E. The problem of predicting sudden events (e.g. catas- 
trophic eruption) then reduces to finding the range, if 
any, in which the slowly evolving structure may suffer 
a loss of equilibrium. The signature of the loss of equi- 
librium is quite easy to identify for the DB states. The 
transition may occur in one of the following two ways: 
1) when the roots of the quadratic equation, determining 
the length scales for the field variation, go from being 
real to complex (implying change in the topology of the 
magnetic and the velocity fields — boundary separating 
the paramagnetic from the diamagnetic), or 2) the am- 
plitude of either of the two states ceases to be real. For 
our current problem, the sudden change is likely to follow 
the second route. 

By analysing a simple analytically tractable model, 
we find affirmative answers to all the four questions we 
posed. We show that the invariants /ii,/i2 , and E , 
which label and alongwith the initial and bounbary con- 
ditions determine the original state, hold the key to the 
eventual fate of a structure. If for a given equilibrium 
sequence, the total energy E is larger than some critical 
value (given in terms of hi, and /i2), the catastrophic loss 
of equilibrium could certainly occur. The trigger for the 
equilibrium loss could come, for instance, from nearby 
structures getting close to each other with an increase in 
their interaction energy. The catastrophy pushes a DB 
state to relax to a minimum energy single Beltrami field. 
For coronal structures, the transition transfers almost all 
the short-scale magnetic energy to the fiow energy. 



II. MODEL 

Within the framework of our approach, there are 
two distinct scenerios for eruptive events : a) when 
a "slowly" evolving structure finds itself in a state of 
no equilibrium, and b) when the process of creating a 
long-lived hot structure is prematurally aborted; the 
flow shrinks/distorts the structure which suddenly shines 
and/or releases energy or ejects particles. The latter 
mechanism requires a detailed time-dependent treatment 
and is not the subject matter of this paper. The following 
semi-equilibrium, coUisionless magneto-fluid treatment 
pertains only to the former case. 

A given structure is supposed to corresspond to the 
equilibrium solutions of the two-fiuid system 



d_ 

dt 



- V X {Uj 



(j = l,2) 



(1) 



written in terms of a pair of generalized vorticities 
u)i = B, UJ2 = B + V X V , and effective flows Ui = 
V — W X B, U2 = V, with the following normalizations: 
the magnetic field B to an appropriate measure of the 
magnetic field Bq, the fiuid velocity V to the correspond- 
ing Alfven speed, and the distances to the coUisioUess ion 
skin depth li. 

The simplest and perhaps the most fundamental equi- 
librium solution to (|l|) is given by the "Beltrami condi- 
tions" , which imply the alignment of the vorticities and 
the corresponding flows {ujj/ /Uj), 



B = a{V -S/ X B), 
B + \7 xV = bV, 



(2) 
(3) 



with a, b — const. We have used constant density as- 
sumption for simplicity - extension to varying density is 
straightforward Equations (||) and (||) combine to 
yield: 



(curl - A+)(curl - A_)B = 0, 
where a — l/a , W x — "curl" , and 
1 



^- = 2 



(6 - a) ± v/(& + a)2 -4 



(4) 



(5) 



For sub-alfvenic flows (the flows we generally en- 
counter in the solar atmosphere), the length scales (A^^) 
are quite disparate. We assume A^jl^ >> Al^ without loss 
of generality. The general solution to the "double Bel- 
trami equations" (^ is a linear combination of the single 
Beltrami fields G± satisfying (curl — X)G = 0. Thus, for 
arbitrary constants C±, the sum 

B ^C+G+ + C-G- (6) 

solves (^, and the corresponding fiow is given hy V = 
(A+ + a) C+G+ + (A- + a) C-G-. 
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The DB field encompasses a wide class of steady states 
of mathematical physics - from the force-free param- 
agnetic field to the fully diamagnetic field. The Bel- 
trami conditions also demand "generalized Bernoulli con- 
ditions" which allow pressure confinement when an ap- 
propriate flow is driven ||^ (and references therein). 

The DB field is characterized by four parameters: 
A+, A- (eigenvalues), and C+, C_ (amplitudes). The 
three invariants ||T^: the helicity hi, the generalized he- 
licity /i2, 



= i J{A-B)dr, 



h2 = - I iA + V)-iB + V xV)dr, 



(7) 
(8) 



{A is the vector potential), and the total energy 

E^\j{B^ + V^)dv (9) 

will provide three algebraic relations connecting them 
iQ. To predict the possibility of an eruptive event, 
interpretted as the termination of an equilibrium se- 
quence (for solar flares, this kind of an approach, albeit 
in different contexts, has been followed in numerous in- 
vestigations, (see e.g. and references therein), 
we analytically investigate this system using the macro- 
scale (A^^) of the closed strcuture as a control param- 
eter. This choice is physically sensible and is moti- 
vated by observations because in the process of structure- 
structure interactions, "initial" shapes do undergo defor- 
mations/distortions with rates strongly dependent on the 
initial and boundary conditions. 

For simplicity we explicitly work out the system in a 
Cartesian cube of length L. We take G± to be the simple 
2-D Beltrami ABC field ITtI, 



G± = gx± I sinAiX 
cos X±x 



cos \±y 
9y± I I , (10) 
sin A±y 



with the normalization {gx±Y + {9y±Y — 1- For real A±, 
( p^ ) represents an arcade-magnetic field structure resem- 
bling interacting coronal loops [in Fig. 1] . Assuming L = 
n+(27r/A+) = n_(27r/A_) (n± are integers), G± satisfy 
the following relations: / G±dr — L^, / G+ ■ G dv = 0, 

where / dr = J^^ dxAy. 

The invariants can now be readily evaluated and the 
results can be displayed in several equivalent forms. We 
find the following three equations to be the most con- 
venient for further analysis (/12 — hi + h2, h2 — 
bE- X+X^hi): 



ho = — 



(A+ + A_) ± v/(A+-A_)2+4 

^ L 

— A_|-A_/ii, 

Cl = D-^ {^^ - [1 + (A- + &f] X-hi} A+, 
= -D-^ {^^ - [1 + (A+ + af] X+hi} A_, 



(11) 
(12) 

(13) 



where we have removed the common factor L'^ /2, and 

D=[1 + {X++ af] A+ - [1 + (A_ + 5)2] A_ 
= {X+- X^)b{b + a). (14) 

For given hi, E, /12 (/12) and A+ (control parameter), we 
can solve the preceding system to determine the phys- 
ical quantities A_, and C± which must all remain real 
for an equilibrium. Before we give an analytic derivation 
for the bifurcation conditions (leading to loss of equilib- 
rium), we display in Fig. 2 the plots of A_ and C± as 
functions of A+ for two distinct sets for the values of the 
invariants: we choose hi = 1, /12 = 1-5, E = 0.4 for 
Fig.2(a), and hi = 1, /12 = 1.5, E = 1.3 for Fig.2(b) 
(dashed lines correspond to the region of imaginary C_). 
We find that the behavior of the solution changes dras- 
tically with E. For the parameters of Fig. 2(a), A_ and 
C± remain real and change continuously with varying A+ 
implying that as the macroscopic scale of the structure 
(1/A+) changes continuously, the equilibrium expressed 
by ( |lO|) persists - there is no catastrophic or qualitative 
change. For Fig. 2(b) with E changing from 0.4 to 1.3 
(with same hi, /i2) we arrive at a fundamentally differ- 
ent situation; when A+ exceeds a critical value A"'*, i.e, 
the macro-scale becomes smaller than a critical size, the 
physical determinants of the equilibrium cease to be real; 
the sequence of equilibria is terminated. 

The condition for catasrophe turns out to be a con- 
straint involving hi, /12 and E, which will allow the van- 
ishing of for positive C^, and real A±. It is straight- 
forward to show that the system has a critical point if 



£:2 >El^A [hi ± ^/hih2 



(15) 



and the critical A is determined by a simultaneous solu- 
tion of (11) and E - [l + (A+ + a)^] X+hi = giving: 



Al"* = 



2hi 



E ± _ El 



(16) 



Thus, for E > Ec (determined by helicities hi and /12), 
when the macroscopic size of a structure shrinks below a 
critical value, it can go through a severe reorganization. 

At the critical point, an expected but most remarkable 
transition occurs. Using the value of A"'*, we find from 
equation(|l^) that the coefficient C_ , which measures the 
strength of the short scale fields, identically vanishes, and 
the equilibrium changes from Double Beltrami to a single 
Beltrami state defined by A+ — A^j5'*, i.e., B = C+G+ 
{V X B = X+B) with V parallel to B. The transition 
leads to a magnetically more relaxed state with the mag- 
netic energy reaching its minimum with appropriate gain 
in the flow kinetic energy (see Fig.||). 



III. CONCLUSIONS 

By modelling quasi-equilibrium, slowly evolving coro- 
nal structures as a sequence of Double-Beltrami mag- 
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netofiuid states in which the magnetic and the veloc- 
ity field are self-consistently coupled, we have shown 
the possibility of, and derived the conditions for catas- 
trophic changes leading to a fundamental transformation 
of the initial state. The critical condition comes out as 
an inequality involving three invariants of the coUisionless 
magnetofluid dynamics. When the total energy exceeds 
a critical energy the DB equilibrium suddenly relaxes to 
a single Beltrami state corresponding to the large macro- 
scopic size. All of the short-scale magnetic energy is lost 
having been transformed to the flow energy and partly 
to heat via the viscous dissipation of the flow energy. 

This general mechanism in which the flows (and their 
interactions with the magnetic field) play an essential role 
could certainly help in advancing our understanding of a 
variety of sudden (violent) events in the solar atmosphere 
like the flares, the erupting prominences, and the coro- 
nal mass ejctions. The connection of flows with eruptive 
events is rather direct: it depends on their ability to de- 
form (in specific cases distort) the ambient magnetic field 
lines to temporarily stretch (shrink, destroy) the closed 
field lines so that the flow can escape the local region with 
a considerable increase in kinetic energy in the form of 
heat/bulk motion. 
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FIG. 1. Magnetic field line structure of a 2-D ABC map 
resembling coronal arcades. 

FIG. 2. (a). Plots of A_ and C± versus A+ for 
E = 0.4 < Ec — 0.45, the critical energy. No catastrophe, 
(b). Plots of A- and C± versus A+ for £ = 1.3 > . There 
is a critical point at A+ ~ 0.041. 

FIG. 3. Plots for A_, the magnetic and the flow energies 
versus A+ for the catastrophe-prone set /ii = 1, /12 = 1.007 
and E = 1.?, > Ec = 7 ■ 10"^. The scale lengths are 
highly separated A+ << A_. The initial choice makes 
C+ ~ 0(A+/A_) < 1 and C- ~ 0(A-/A-) ~ 1 from (|l|) 
and (p^. If any interaction increases A+ (the size of the struc- 
ture shrinks) the critical point (A+ = A+'') will be reached at 
which C- is zero. The magnetic field energy (oc C\ + C?.) 
decreases to a very small value since Cl < 1. Since the 
total energy is conserved, almost all the initial magnetic en- 
ergy is transfered to the fiow causing an eruption. Notice 
that for coronal plasma, the skin depth U is small ~ 100cm 
{h/\+ ~ 10"^ km), for a typical density of ~ 10^ cm~^. A word 
of caution is necessary - as we approach the critical point, the 
quasiequilibrium considerations are just an indicator of what 
is happening - they must be replaced by a full time-dependent 
treatment to capture the dynamics; the changes are no longer 
slow. 
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